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Abstract. The astrophysical S(E) factor of 14 N(p,7) ls O has been measured for effective center-of-mass 
energies between E e //= 119 and 367 keV at the LUNA facility using TiN solid targets and Ge detectors. 
The data are in good agreement with previous and recent work at overlapping energies. R-matrix analysis 
reveals that due to the complex level structure of 15 O the extrapolated S(0) value is model dependent and 
calls for additional experimental efforts to reduce the present uncertainty in S(0) to a level of a few percent 
as required by astrophysical calculations. 

PACS. 25.40.Lw Radioactive capture - 26.20.+FHydrostatic stellar nucleosynthesis - 26.20.+tSolar neu- 
trino - 29.30.KvX - and 7 ray spectroscopy - 97.10.CvStellar structure and evolution 



1 Introduction 

The capture reaction 14 N(p,7) 15 0, the slowest process in 
the hydrogen burning CNO cycle PQ, is of high astrophys- 
ical interest as its reaction rate influences sensitively the 
age determination of globular clusters [2] and the solar 
neutrino spectrum [30. The capture cross section needs 
to be known down to Eq = 30 keV (the Gamow peak in 
core H-burning stars), which is far below the low-energy 
limit of direct 7-ray measurements, i.e. the center-of-mass 
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energy E = 240 keV a 0. Thus, the data had to be ex- 
trapolated over a large energy gap leading to a substantial 
uncertainty for the astrophysical S-factor at zero energy, 
S(0). According to the data and analyses of Schroder et 
al [5], there are two major and nearly equal contributions 
to S(0): the direct capture (DC) to the 6.79 MeV state in 
15 and the capture to the ground state (gs) in ls O. The 
latter process is enhanced due to a subthreshold resonance 
at Ex = —507 keV (Fig.QJ, the width of which was taken 
as a free parameter in the fit [2] • The extrapolation led to 
5 tot (0) = 3.20 ± 0.54 keVb, with S gs (0) = 1.55 keVb and 
#6.79(0) = 1.41 keVb. 

Subsequently, the data of Schroder et al [5] were ana- 
lyzed by Angulo and Descouvemont 6 using an R-matrix 
approach [7] . Contrary to the extrapolation by Schroder et 
al [HJ for capture to the ground state, they reported a neg- 

a In this work all proton energies are taken in the center-of- 
mass system, except where quoted differently. 
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ligible contribution S gs (Q) = 0.08 keVb due to a smaller 
total width of the subthreshold resonance and thus they 
suggested 5 tot (0) = 1.77 ± 0.20 keVb. A smaller width of 
the 6.79 MeV state was supported by a lifetime measure- 
ment via the Doppler-shift method 8 and by a Coulomb 
excitation measurement [0]. 

Both S gs (0) and 56.79(0) are dominated by El capture 
mechanisms. Based on a measurement of the analysing 
power at E p = 270 keV, Nelson et al. JU] reported a small 
contribution (about 4%) of an Ml capture to the overall 
capture process, but this single data point did not signif- 
icantly improve the situation. The available data clearly 
demonstrate the large uncertainty affecting the extrapo- 
lations over a large energy gap in a nucleus, such as 15 0, 
with a complex level structure, a problem calling for new 
direct 7-ray measurements towards lower energies. 

The LUNA (Laboratory Underground for Nuclear As- 
trophysics) collaboration started in 2001 a reinvestigation 
of 14 N(p,7) 15 (using TiN solid targets and Ge detec- 
tors) with a particular emphasis on S gs (E). The data ob- 
tained down to E p = 140 keV JT] confirmed the small 
contribution of the ground state capture giving Stot (0) = 
1.7 ± 0.1 ± 0.2 keVb. The R-matrix analysis of these low- 
energy data included the high-energy data of Schroder et 
al. [5] and the results of the ANC (Asymptotic Normal- 
ization Coefficient) method |12H13| . At the completion of 
the LUNA work, another set of low-energy data became 
available at E p — 155 to 524 keV which is in excellent 
agreement with the LUNA data in the overlapping energy 
range. Here, we report on the LUNA solid target work, 
while additional details can be found in |14II15II16| . 



2 Experimental apparatus 

2.1 Accelerator 

The 400 kV LUNA accelerator [TH] provided a proton 
beam on target of up to 500 /xA. The absolute beam energy 
was known with an accuracy of 300 eV, and the energy 
spread and the long-term energy stability were observed 
to be 100 eV and 5 eV/h, respectively. The beam current 
on target was monitored by a current integrator with an 
estimated uncertainty of 2 %. The targets were directly 
water-cooled. In the target chamber (Fig. [2J) the target 
ladder was surrounded by an electrically insulated colli- 
mator biased to -300 V to suppress secondary electrons. 
A uniform rectangular beam spot (4 by 4 cm) was pro- 
duced within the target area by magnetic wobbling of the 
beam. In order to prevent build-up of impurities on the 
target, a LN 2 -cooled Cu cold finger was used. 

2.2 Detectors and electronics 

The laboratory's 1400 m rock cover reduces the back- 
ground rate in the Ge detectors by more than three orders 
of magnitude at E~ > 5 MeV 19 compared with a de- 
tector placed at the surface of the earth. It is in this E 1 
region that the capture 7-ray lines of secondary transitions 



(5.18, 6.17, and 6.79 MeV) and the ground state transition 
(about 7.5 MeV) are located (Fig.QJ. In the experiments 
two different setups were used. 

Geometry 1 (Fig. [21 : The target was oriented with its 
normal at 55° with respect to the beam direction. The 
capture 7-rays were observed with a HPGe-detector 
(126% relative efficiency), whose front face and target 
were parallel with a distance d = 1.5 cm. The detec- 
tor was surrounded by 5 cm of lead, which reduced 
the background in the low-energy range of the primary 
transitions by a factor of 10. The setup was used for 
the measurement of excitation functions, while for the 
determination of the detector efficiency and of sum- 
ming effects the distance between target and detector 
was increased in discrete steps up to d = 20.5 cm. 

Geometry 2: The target normal was oriented at 19° 
with respect to the beam direction. Using a different 
target chamber three HPGe-detectors were positioned 
around the target: the 126% detector was placed at 0°, 
while the other two detectors (120% and 108% relative 
efficiencies) were placed at 90° and 125°. The distance 
between the detector end cap and the target was 7 cm 
for all detectors. The setup was used in the measure- 
ment of Doppler shifts, excitation energies, and angu- 
lar distributions. 

2.3 Targets 

The 14 N solid targets had to fulfill several requirements: 
(i) uniform depth distribution starting at the surface, (ii) 
thickness large enough that possible sputtering effects at 
low energies and over long running times were negligi- 
ble for the yield measurements, (iii) stability against high 
beam power, and (iv) low concentration of impurities of 
light elements. The targets were produced either by im- 
plantation of a 14 iV ion beam into a suitable backing, or by 
evaporation of a Ti layer on a specific backing and heat- 
ing it in a N2 atmosphere, or by sputtering of TiN on a 
backing. The resulting quality of the targets was investi- 
gated using the Er — 259 keV resonance. Fig.|3]shows the 
results out of a large number of tests |14j performed with 
different backings and target thicknesses. The evaporated 
and sputtered targets showed a steep rise in the excitation 
function at the resonance energy, with a slope compara- 
ble to the convolution of beam resolution and resonance 
width. In the implanted target the surface region was not 
saturated. The step heights to the plateau were similar 
for all targets indicating a similar stoichiometry at satura- 
tion. Our conclusion was that the sputtered targets of TiN 
(with a typical thickness of 80 keV, Ta-backing) had the 
most uniform number density profile and could withstand 
many days of intense beam bombardment without a sig- 
nificant deterioration. Typically, a new TiN target was in- 
stalled after a running time of 1 week, which corresponded 
to an accumulated charge of about 200 C. During this time 
we observed a small decrease of the plateau-yields due to 
the proton implantation (insert of Fig. [3| • The stoichiom- 
etry of the TiN layer - used for the determination of the 
strength of the Er — 259 keV resonance - was measured 
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Fig. 1. Relevant level scheme of 15 O near the 14 N(p,7) 15 Q threshold. 
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Fig. 2. Schematic diagram of the setup used in geometry 1. 
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via Rutherford Backscattering Spectrometry using a 2.0 
MeV 4 He beam in the INFN National Laboratory of Leg- 
naro, resulting in Ti/N — 1/(1.08±0.05). To measure the 
capture cross section above the Er = 259 keV resonance 
we used evaporated targets, since we needed here thinner 
targets to minimize the influence of the resonance: their 
absolute stoichiometry did not need to be precisely known 
because of the normalization procedure adopted (section 

E3J. 



2.4 Beam induced background 

The presence of light-element impurities can produce in- 
tense 7-ray background lines due to their relatively low 
Coulomb barrier and/or high cross sections. In particular 
7-lines from reactions on 11 B, 12 C, 13 C, ls O and 19 F have 
been investigated and their intensity was minimized by a 
proper choice of target backing and preparation proce- 
dures where the sputtered target was best. A sample 
7-ray spectrum obtained at E p = 250 keV is shown in 
Fig. 0] the primary capture transitions (e.g. tr — * 6.79) 
and the secondary transitions (e.g. 6.79 — > 0) as well as 
background lines (identified by their impurity nuclides) 
are indicated. The high-energy part of the spectrum ob- 
tained at E p = 140 keV, our lowest bombarding energy, is 
shown in Fig. [S] the primary transition tr — > and nearly 
all secondary transitions are clearly observable. 

3 Experimental procedures, data analyses and 
results 

3.1 Excitation energies in 15 

The energy scale of the Ge detectors used in geometry 
2 was determined using the well-known background lines 
from 232 Th, 40 K and 214 Bi and the ground state capture 
transition at the resonance. The latter energy was calcu- 
lated from the known accelerator energy (Er = 259.4±0.3 
keV, section ESI and the Q value of 14 N(p, 7 ) 15 (Q = 
7296.8 ± 0.5 keV J2H), where effects of the Doppler shift 
and recoil energy were taken into account: E x = 7556.4 ± 
0.6 keV. The energies of the 77 cascades via the 5.18, 
6.17 and 6.79 MeV states were then used to determine 
the excitation energies of these states, where the Doppler 
attenuation factors were left as a free parameter in the 
fit. The results are summarised in tableland are in good 
agreement with previous work [21| . 



3.2 Efficiency and summing effects 

Due to (i) the close geometry of the set-up in geometry 
1 (Fig. 2), (ii) the cascade structure of the capture pro- 
cess, and (hi) the small branching to the ground state, 
the "summing-in" effect (i.e. the cascade 7-rays giving a 
contribution to the ground state full-energy peak) and the 
" summing-out" effect (i.e. a full-energy detection of one 



cascade 7-ray and any concurrent interaction of the other 
member of the cascade) have to be taken into account. For 
the efficiency and summing-effect evaluation, we consid- 
ered 7-ray spectra obtained at 55° for different distances 
d from the target (i.e. d = 1.5, 5.5, 10.5, and 20.5 cm). The 
spectra involved 137 Cs and 60 Co calibrated sources placed 
at the target position as well as the spectra obtained at 
the Er — 259 keV resonance. The dependence of the full- 
energy efficiency Sf e on the 7-ray energy E~ l and on the 
distance d were parameterized by the following functions 

ln(e /e ) =a + 61n(£; 7 )+c[ln(£; 7 )] 2 (1) 

d + d.Q 

£fM = (d + d ) 2 ® 

In order to take into account the summing-in and summing- 
out effects, the following expressions were used 

Y gs = R(b gs e fe (E gs )+Y / b l e fe (Er)efe(En) 

Y ipri = Rbi£ fe (E ipri )(l - etot(E ls J) 

Yisec = Rh£ f e(E ls J(l - e tot (E ipri )) (3) 

where Y gs , Yi pri and Yi Bec are the observed yields of the 
ground state transition and of the transitions through the 
energy level i; the subscripts pri and sec refer to primary 
and secondary transitions, respectively; R is the number 
of reactions per unit charge, h gs and hi are the branchings, 
£f e and Etot are the full energy and total efficiencies, re- 
spectively. The equation 01 assumes isotropic angular dis- 
tributions and correlations. In these equations the total 
detection efficiency is parameterized as |23| : 

l n (££l) = K 1 +K 2 ln(£ 7 ) + if 3 (ln(£g) 2 (4) 

Etot 

Each excited state that is fed by the resonance then decays 
to the ground state with a 100% probability: this allows to 
set constraints on the efficiency curve arising from the one- 
by-one equality of the intensities of primary and secondary 
transitions involving each excited state: 

Y lp „ = e fe (E lprt )(l-e tot (E ls J) 

Y iaec e }e (E iaec )(l-e tot (E ipri )) K ' 

These constraints were used in a global fit to the data. 
The free parameters are the coefficients in the parame- 
terizations (equations ^ HI and 0} and the four branching 
ratios of the resonance. Fig. HJ shows the behavior of the 
efficiency £f e (E) versus energy for the four distances: the 
data points are those from the reaction (diamonds) and 
from the sources (squares) and the lines are the fit. Points 
and data refer to distances of 1.5, 5.5, 15.5 and 20.5 cm 
starting from the highest curve, respectively. In order to 
illustrate the importance of the summing correction, the 
crosses in the figure represent the values of the efficiency 
at the closest distance which would be obtained neglect- 
ing the summing effects: note the large deviation for the 
ground state capture at 7.5 MeV. The results shown in 
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Fig. 3. Excitation function of 14 N(p,7) 15 near the Er = 259 keV resonance using different targets: squares = sputtered 
target, triangles = implanted target, diamonds = evaporated target. The insert shows the thick-target yield as a function of 
accumulated charge Q on the sputtered target. 
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Fig. 4. Spectrum obtained at E p = 250 keV in geometry 1 over an accumulated charge of 87 C. The background above U 7 = 8 
MeV is mainly due to a small 11 B contamination. 
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Fig. 5. Spectrum obtained at E p = 140 keV in geometry 1 over an accumulated charge of 210 C. 



Fig. [f)| are consistent with values obtained from Monte- 
Carlo simulations |15| . 

The resulting branchings of the E R — 259 keV reso- 
nance are in good agreement with previous work |2I] (ta- 
bic nj except for the ground state decay. 



3.3 Strength, energy and width of the E R — 259 keV 
resonance 

The thick-target yield Y°° of the resonance is given by the 
expression p^: 

A 2 LU^f 



energy E is given by the integral: 



2 V (E) 



(6) 



where r)(E) = dE/dx is the stopping power of the target 
at the resonance energy. For the observed stoichiometry 
of the TiN target, the compilation [21] leads to a stop- 
ping power t](Er) = 32.8 eV / (10 15 atoms / cm 2 ) . We find 
then a resonance strength uj^ — 12.9 ± 0.4 (stat) ± 0.8 
(sys) meV. The systematic error arises mainly from the 
uncertainty in the stopping power and the accuracy of 
the current measurement. The strength value is in good 
agreement with previous work: ujj = 14 ± 1 meV |21j 
and u>7 = 13-5 ± 1.2 meV Q7]. 

In the energy region of the En — 259 keV resonance, 
the yield of a capture transition as a function of beam 



Y(E) 



dE 



(7) 



With the known stopping power and assuming the cross 
section as an incoherent sum of a constant non-resonant 
term plus a Breit-Wigner term, we have fitted the ex- 
perimental curves of all four primary capture transitions, 
where the free parameters were E^ and r R . The results for 
all four transitions (Fig.0) were within their respective er- 
rors leading to E R = 259.4 ±0.3 keV and r R = 0.96 ±0.05 
keV, in good agreement with previous work |51l21| . 



3.4 Angular distributions 

Angular distributions were measured using geometry 2 at 
energies E = 206, 313, and 361 keV. Since the r = 1/2+, 
E R = 259 keV resonance exhibits isotropic distributions, 
the relative efficiencies of the Ge detectors were measured 
at this resonance. The resulting distributions were fitted 
with Legendre polynomials 



W(6) = 1 + aiQiPi(0) + a 2 Q 2 P 2 (6) 



(8) 



where are the angular distribution coefficients and Qk 
the attenuation coefficients calculated according to |22| . 
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Fig. 6. Full-energy peak efficiency as function of 7-ray energy for the geometry f with distances d — 1 .5, 5.5, f5.5, and 20.5 
cm, from top to bottom. The lines through the data points are the results from a fit. The crosses are the results omitting the 
summing effects. 
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Fig. 7. Excitation functions of the 4 primary transitions near the Er = 259 keV resonance. The lines through the data points 
are the result of a fit. 
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Within experimental uncertainties, all primary and sec- 
ondary transitions are isotropic |15| and consistent with 
the results of Schroder et al. [S]- As known from the 
exception is the primary transition to the 6.79 MeV state, 
which exhibits a sin 2 (#) distribution (a 2 = —0.50): a 2 = 
-0.79 ±0.11, -0.64 ± 0.06, and -0.57 ± 0.07 at E = 206, 
313, and 361 keV, respectively, consistent with the data 
from 5] at higher energies. 

3.5 Primary peaks 

Excitation functions have been obtained at proton ener- 
gies between E p = 140 and 400 keV for all primary tran- 
sitions. The observed line-shape of a primary transition is 
determined by the cross section behavior <r(E p ) in the pro- 
ton energy interval spanned by the incident beam during 
the slowing-down process in the target, once the transfor- 
mation from the energy E p (at which the reaction takes 
place) to the corresponding 7-ray energy E 7 = E+Q-E^ 
is performed. The shape is also influenced by the energy 
loss of the protons in the thick target, since the stopping 
power of the protons in TiN is a function of proton en- 
ergy [25 . Finally, the position of the high-energy rise in 
the line-shape is influenced by the Doppler effect and the 
recoil of the compound nucleus. 

The number of counts Nj in channel i of the 7-spectrum, 
corresponding to the energy bin E 7 i to E 7 i + (5E 7 (<5E = 
dispersion in units of keV per channel) is given by the 
expression 



Ni 



a(E p i)8E 1 ef e (E 1 {)b : j 
ri(E p i) 



(9) 



for E p j < Ep (E p i = proton energy corresponding to chan- 
nel i, E p = incident proton energy), where a(E P i) is the 
capture cross section, Sf e {E li ) is the 7-ray detection ef- 
ficiency, and bj is the branching of the associated decay. 
The conversion from E 7 ; to E p i includes the Doppler and 
recoil effects. The result is folded with the known detec- 
tor resolution AE~ t to obtain the experimental line-shape. 
The spectrum obtained at an energy E p was compared to 
that obtained with the same target at the Er = 259 keV 
resonance energy: the comparison was performed before 
and after a given run at E p , where the thick-target res- 
onance yields Y°° were averaged. In this comparison the 
experimental fit-quantity is the ratio 



2a(E pi )5E 1 e fe {E li )bjri{E R i ab 



using ai(E p ) deduced from eauation llOl 

The summing-in correction for the ground state tran- 
sition is not negligible and becomes increasingly uncertain 
with decreasing proton energy. This was the main reason 
why we stopped the measurements at E p = 140 keV al- 
though from a statistical point of view and background 
considerations we could go lower in energy (Fig. |SJ. In 
the case of the other capture transitions the summing-out 
correction depends only on the precision of the efficiency 
curve. 

The last step in the analysis of the primary peaks is to 
introduce the function for the cross section in eauation llOl 
Since we have analysed the cross section above and below 
the Er = 259 keV resonance, we have written - in first 
approximation - the cross section as an incoherent sum 
of a non-resonant term, which assumed a constant astro- 
physical S(E) factor Snr, and a resonant term described 
by the Breit-Wigner formula: 



\ 2 u!-ybRjefe(E R7 )ri(E P i) ^ ^ 

The energy dependence of the full-energy efficiency is known 
the summing-out correction for the cascade transitions via 
the three excited states is smaller than 1%; the summing- 
in correction for the capture to the ground state depends 
on bj and b R j. Thus, for the ground state capture a cor- 
rection has to be applied which can be obtained replacing 
equation with: 

a gs (E p )s fe (E!>°) + £. ai(E p )s fe (E°?)e fe (E0 



yoo 



oc 



o{E p 



—2irr) 



E,, 



A 2 

7T ^(Ep, 



r 



Er) 2 + (r/2) 2 ) 



* gs (E Rp )e fe (E° R ) + Ei^(^)e/e(^)e/e(£?S 

(11) 



These expressions were fitted to the 7-ray spectrum yields, 
superimposed on an exponential background. 

An example of this fitting procedure is given in Fig. 
03 where the 7-spectra in the region of the ground state 
transition (top figure) and of the primary transition to the 
6.79 MeV state (bottom figure) are shown as obtained at 
Ep = 250 keV: the gray area is the background, the left- 
oriented stripes area and the right-oriented stripes area 
represent the non-resonant and resonant parts of the cross 
section, respectively, and the squared area is their sum. 
In the fits we used as free parameters the non resonant 
astrophysical factor and the background parameters. After 
the R-matrix analysis (section 0} this procedure has been 
checked resulting in a fit shown at the bottom of Fig. 
for the transition to the 6.79 MeV state. 

Since the non-resonant astrophysical S(E) factor can 
be considered as constant only in a short energy range, 
our thick target may represent a too large energy win- 
dow. Indeed the value, where the fitting value converges, 
varies with energy. Moreover, in some cases background 
peaks can be present in the line-shape region, due either 
to natural background or beam-induced background. For 
example, in the region of the ground state peak there are 
contamination peaks due to the presence of 18 O in the tar- 
get, and in the region of the low-energy primaries there are 
many lines due to the natural radioactivity of the rocks. 
Thus, an alternative way to analyse the data, i.e. simu- 
lating a thin target, is to subtract two spectra acquired 
at different energies. The resulting spectrum is equivalent 
to one acquired at the higher energy with a target thick- 
ness equal to the difference in incident energy (A = E max 
- E m i n ). An example is shown in Fig. [5] the upper part 
shows the 7-spectra obtained at E p = 250 keV (dark gray 
area) and 240 keV (light gray area) for the primary tran- 
sition to the 6.79 MeV state; the center part shows the 
-7-spectrum resulting from the subtraction of both spec- 
)tra, which is equivalent to a spectrum acquired at an en- 
ergy of 250 keV with a target of A = 10 keV thickness. 
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Fig. 8. Gamma-ray spectra near the ground state transition (top) and the primary transition to the 6.79 MeV state obtained 
at E p = 250 keV (for details, see text). 



After the final R-matrix analysis was completed we have 
reinvestigated this differential method, using the R-matrix 
results, which is shown in the bottom part of the Fig|5| 
Excellent agreement with the approximation shown in the 
upper and center parts is noted. The differential method 
in combination with the line-shape analysis has been used 
for the primary transitions. 



3.6 Secondary transitions 

The analyses of the secondary transitions have several ad- 
vantages: (i) since the 7-energies of all secondary transi- 
tions are above 5 MeV, we can make full use of the ad- 
vantage of an underground laboratory; (ii) the expected 
peaks are relatively narrow (compared to the primaries) 
and hence the peak to background ratio is improved; (iii) 
there is no efficiency correction necessary when comparing 
the off-resonance yield with that on-resonance; (iv) the in- 
terpretation of the observed yield in combination with an 
80 keV thick target is nearly identical with the integral 
of the cross section from E=0 to the incident energy, thus 
target-profile corrections are not necessary; (v) angular 



distribution effects are expected to be negligible (section 

roii . 

The yields of all secondary transitions were obtained 
by fitting the background with a linear function at the low- 
energy and high-energy sides of the peak which was then 
subtracted from the total peak intensity. The yield was 
then compared with the corresponding thick-target yield 
step Y°° at the resonance. The resulting thick-target yield 
curves are shown in the inserts of Figs. II f 1 1121 and ITU for 
the secondary transitions E 7 = 5.18, 6.17, and 6.79 MeV, 
respectively. Each yield ratio corresponds to the integral 
from E = to the respective incident beam energy E< 
taken in the center of mass system. EauationllOlcan thus 
be written: 



a{E)b P1 {E R ) 



dE 



(13) 



where uj^f are the respective values for the individual 
transitions of the E^j = 259 keV resonance. These yield 
ratios are corrected for summing-out effects. Note that the 
stopping power only enters equation^] as a ratio and has 
not to be known absolutely. In principle one could also 
use the differential method (subtraction of the secondary 
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peak contents from two runs with successive beam ener- 
gies) to determine the energy behavior of the cross section. 
However, the statistical error increases by subtracting two 
large numbers together with the uncertainty in the deter- 
mination of the effective beam energy. Since the energy 
dependence of the cross section is in principle known from 
the R-matrix analysis of the primaries we have adopted 
the method of fitting the thick target yield curve (equa- 
tion ii3i> using the R-matrix analysis with a fine tuning 
of the fit parameters. The solid curves through the data 
points of the secondary transitions, the insert in Figs. lTTl 
I12land fm are the results of integrated R-matrix calcula- 
tion. 

To complete the study of the secondary peaks we have 
studied the 5.24 MeV (see fig.[T3J)transition. The 5.24 MeV 
state is populated by cascades through the 7.28 and 6.85 
MeV states, which in turn decay 100 % to the 5.24 MeV 
state 5 . In the 259 keV resonance the 5.24 MeV state 
is populated directly with 0.6 % branching. At energies 
below this energy the primary peaks feeding the 7.28 and 
6.85 MeV states were not visible in our spectra. 

The results of both analyses are shown in figures ITOI 141 
and numerical values are given in tabled The solid curves 
in the insert of figures 1101141 are the results of integrated 
R-matrix calculations. 



4 R-matrix analysis and discussion 

The experimental results for the transition to the ground 
state and to the = 6.79 MeV state were reported in 
[TT] and an R-matrix [7||nj analysis was described for these 
two transitions with the results: S ffs (0) = 0.25 ±0.06 keVb 
and S^.rg = 1.35 ± 0.05 keVb, where the quoted errors in- 
clude only statistical uncertainties of the data; this error 
quotation was also applied below. The cross section for all 
transitions to excited states can be determined via the pri- 
mary and secondary transitions. The S(E) determination 
of the primary transitions has been described in section 
13. 51 and can be directly used for an R-matrix analysis. The 
peak of a secondary transition for an infinitely thick tar- 
get contains the integrated yield from zero energy to the 
respective beam energy. Thus, its interpretation requires 
knowledge of the energy dependence of the cross section. 
This energy dependence is provided by the R-matrix fit 
on the primary data as a starting condition with a sub- 
sequent fine tuning of the fitting parameters. Such a pro- 
cedure was used for the transitions to the 5.18, 6.17 and 
6.79 MeV states. The best choice of the channel radius for 
fitting the transition to the ground state and 6.79 MeV 
state was found to be a = 5.5 fm, which was used for 
the present R-matrix calculations. 

4.1 Transition to the ground state 

The data from LUNA and previous works for the ground 
state capture are shown in FigE3 We refer to the analysis 
described in our previous paper [TT] . Briefly, the data were 
fitted including the 3/2 + sub-threshold state, the 1/2+ , 



259 keV, the 3/2+, 987 keV and the 3/2+, 2187 keV res- 
onances as well as a background pole located at 6 MeV. 
For the sub-threshold state, we used the reduced width 
obtained from the fit of the data for the 6.79 MeV transi- 
tion. The fit parameters were the TV of the subthreshold 
state, the 987 keV and 2187 keV resonances, the r p and 
iV of the background pole. For the external contribution 
we used the ANC of ^3] as a starting value. The corre- 
sponding r~f value for the subthreshold state is 0.8 ± 0.4 
eV, which is in good agreement with the value from a life 
time measurement by 8 T 7 = OTI+q'^ eV as well as with 
T 7 = 0.95+^°, the value fr om coulomb excitation work 
[Sj. The minimum in the S-factor shape at energy lower 
than 200 keV is due to an interference effect between the 
sub-threshold state and the nonresonant capture mech- 
anism. The extrapolation of the S-factor is based on the 
energy dependence of the complete data set and not on the 
lowest data points. The extrapolated value of the ground 
state S-factor is S gs (0) = 0.25 ±0.06 keVb. Since the pub- 
lication of JJJ a new experiment was reported ^3 with 
an extrapolated S gs (0) = 0.49 ± 0.08 keVb (table 0). The 
present data and those of are in excellent agreement 
within their respective errors but the extrapolated S 9S (0) 
differ by a factor of two however still within their 2a errors 
(Fig. Il()|l . The major difference in both R-matrix analyses 
is that in ^JJ the high energy data of [5] - corrected for 
summing-in effects - were included in our analysis forcing 
the Sg S (E) to be higher above the 259 keV resonance and 
lower below the resonance. 



4.2 Transition to the E K =6.79 MeV state 



It was found that the % 2 of the R-matrix fit for the tran- 
sition to E x — 6.79 MeV state was primarily determined 
by the data of |5] forcing the fit to be higher than the 
lowest energy data points of JJJ. It was thus desirable to 
check the data above the En = 259 keV resonance. This 
experiment has been performed recently at the 4 MV Dy- 
namitron tandem accelerator in Bochum over the energy 
range E p = 0.6 to 1.3 MeV \fi^ using the setup of geome- 
try 1 (Fig. EJ) • The results are used together with the data 
of jll| in a new R-matrix fit and are shown in fig. 1111 The 
high-energy data are lower than those of [5] but still within 
the respective 2a errors. In this case the fit parameters for 
the R-matrix analysis were 1^ of the En — 259 keV reso- 
nance and the ANC of the 6.79 MeV state. The resulting 
fit of the thick target yield curve for the E 1 = 6.79 MeV 
secondary peak is also shown as insert in fig llll The new 
extrapolated value S6.7g(0) = 1.21 ± 0.05 keVb is within 
2a of the value given in 1 1 1 1 and in excellent agreement 
with S 6 . 79 (0) = 1.15 ±0.05 keVb from Q7|. The new ANC 
value of C — 4.7 ± 0.1 fm -1 / 2 is still in excellent agree- 
ment with the ANC results Q3| (C = 5.2 ± 0.7 fm" 1 / 2 ) 
and U (C* = 4.6 ± 0.5 fm" 1 / 2 ). Also shown in figdl is 
the R-matrix fit of JOj as a dashed curve: the fit does not 
reproduce the present data. 
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Fig. 10. The astrophysical S 9S (E) factor for the capture to the ground state is shown from the LUNA studies (filled- in points), 
from previous work [Sj (open squares) and from the LENA group )17) (crosses). The thick line through the data points represents 
our R- matrix fit, while the dashed line represents the LENA R- matrix fit. 



4.3 Transition to the Ea;=5.18 MeV state 

This transition can be fitted with the inclusion of the two 
,r = 1/2+ resonances at E R = 259 and 1446 keV and the 
J w = 3/2+ resonance at E# = 987 keV, which is added 
incoherently. The result is shown in fig.^]for the primary 
and secondary transitions. A small external DC contribu- 
tion could be considered but has been neglected here since 
it amounts to less than 1% to the total S-factor. The re- 
sulting S s .i 8 (0) = 0.010 ± 0.003 keVb can be compared 
with 0.018 ±0.003 from 0. 

4.4 Transition to the E.^5.24 MeV state 

The E 7 = 5.24 MeV secondary transition was observed 
with a similar low intensity as the E 7 = 5.18 MeV tran- 
sition. Due to the lack of detailed data, a R-matrix fit 
was not attempted and only a fit of an exponential rep- 
resenting the external DC contribution was used (Fig. 
H3J) yielding S 5 . 2 4(0) = 0.070 ± 0.003 keVb compared to 
S 5 . 24 (0) = 0.064 ± 0.002 keVb from 0. 

4.5 Transition to the E.^6.17 MeV state 

The fitting of the transition to the 6.17 MeV state (Fig. 
114(1 requires an external DC contribution together with 



the two J 71 " = 1/2+ resonances at E R = 259 and 1446 
keV which all interfere with each other, together with the 
J T = 3/2+ resonance at Er = 987 keV added incoherently 
(Fig. ^) . Again, the data above E p = 350 keV are from 
[S] . The best fit was obtained with the respective resonant 
r 7 widths 22.6, 70, and 10 meV. In particular, the fit at 
energies between the two J v = 1/2+ resonances is not 
very sensitive to the choice of the ANC but it is sensitive 
to this choice at the low-energy wing of the E/j ; = 259 keV 
resonance. The yield of the secondary transition (insert 
in Fig. HU) is fitted best with C = 0.2 ± 0.1 fin" 1 / 2 but 
is not acceptable if C exceeds 0.3 fm -1 / 2 . The result is 
within two standard deviation of .13; C — 0.47 ± 0.03 
f m -i/2 anc j f r^| (j — 0.45 ±0.05 fm -1 / 2 after conversion 
to the present coupling scheme. The extrapolated values 
S 6 .i 7 (0) = 0.14 ± 0.02 5 and 0.13 ± 0.02 keVb [Ej are 
consistent with the present result 56.17(0) = 0.08 ± 0.03 
keV b. Recently, the analysing power was measured jl()| 
at E p = 270 keV; the result together with the data of [H] 
led to Sq. 17(0) = 0.16 keV b, which is a factor two higher 
than our value. The agreement of their fit (scanned from 
|10j) with our data at the low-energy wing of the 259 keV 
resonance is quite good but seems to become too high at 
the high-energy side of the 259 keV resonance where the 
influence of the Ml component increases. Unfortunately 
the authors ^Uj only show their fit up to E p = 330 keV so 
that further conclusions are not possible. The analysis of 
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Fig. 11. The astrophysical Se. 79(E) factor for the capture to the 6.79 MeV state is shown from the LUNA studies (filled-in 
points), recent work |16| (open squares)and from the LENA group |17| (crosses). The insert shows the thick-target yield data 
for the 6.79 MeV secondary transition. The lines through the data points represent the results of our R-matrix fit. The dashed 
curve is scanned from |10| . 



[HI with S 6 .i 7 (0) = 0.04 ± 0.01 keVb is a factor two lower 
than the present result. Our data extend towards lower 
energies than those of |17j and in our R-matrix analysis 
an external DC contribution was included which interferes 
with both J v = 1/2+ resonances; thus, no background 
pole needs to be considered. Taking all extrapolations just 
discussed into account, the spread in 86.17(0) is of the 
order of ± 0.06 keVb which corresponds to 4% of the total 
S(0)-factor. 



4.6 Total S-factor 

Our present total S-factor, based primarily on R-matrix 
fits yields S tot (0) = 1.61 ± 0.08 keVb. This result is lower 
by 6% than the value given in fJJ] mainly due to the re- 
vised analysis of the transition to the 6.79 MeV state. The 
present result is in good agreement with |13II17| but dif- 
fers in the weight of the contributions from the various 
transitions (table |3J). 



4.7 Reaction rates 

We obtained an improved extrapolation for S{ O t(0) (table 
121 on the basis of the present R-matrix analysis where all 



presently available experimental information were consid- 
ered within their two a limits. The total reaction rate for 
the 14 N(p, 7) 15 reaction was calculated numerically from 
the R-matrix results with |26| : 

N A < av > = 3.7313 • 10 V*?^ • 

/ S(E)e(- 2 ^e( um5E ^dE (14) 
Jo 

in units of cm^mole -1 s _1 , where S(E) is in MeVb and E 
in MeV and [i is the reduced mass. The resulting rates 
can be expressed within 0.5% using the following analytic 
expression for Tg < 2: 

N A <av>= ai 10 7 T~ 2/3 exp [(a 2 T g - 1/3 - (T 9 /a 3 ) 2 )] • 
(a 4 + a 5 T 9 + a 6 T g 2 + a 7 T 9 3 ) + 
aglO^^expIagTg- 1 ] + 
aiolO^g^exp^Tg- 1 ] (15) 

where: 



ax 


= 3.12 


a 5 


= -1.50 


a 9 = 


= -2.998 


a 2 


= -15.193 


(7,5 


= 17.97 


OlO 


= 8.42 


«3 


= 0.486 


a 7 


= -3.32 


an 


= 0.0682 


04 


= 0.782 
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= -4.891 
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Fig. 12. The astrophysical Ss.i8(E) factor for the capture to the 5.18 MeV state is shown from the LUNA studies (filled-in 
points) and previous work |5] (open squares). The insert shows the thick-target yield data for the 5.18 MeV secondary transition. 
The lines through the data points represent the results of our R-matrix fit. 



a i 


= 2.76 


a 5 


= -1.40 


a 9 


= -2.998 


r .low 
a 2 


= -15.193 


a 6 


= 15.82 


r ,low 
a 10 " 


= 8.44 


~low 
a 3 


= 0.503 


-low 


= -3.32 


a ll - 


= 0.0682 


„low 
u<4 


= 0.804 


~low 

a 8 


= 2.03 


~low 
a 12 ~ 


= -4.987 




high 

a 1 


= 3.44 


h'lqh 
a 5 


= -1.59 


hiqh 

a 9 


= -2.997 


high 
a 2 


= -15.193 


hiqh 

a 6 


= 19.83 


hiqh 

a w 


= 8.42 


high 
a 3 


= 0.475 


hiqh 


= -3.30 


hiqh 

a u 


= 0.0681 


high 
0-4 


= 0.771 


hiqh 

a 8 


= 2.18 


high 
a 12 


= -4.807 



The parameters labeled "low" and "high" in equation 
H5l correspond to the two a limits of present S tot (0) extrap- 
olation and include the error in the strength determination 
W7 of the 259 keV resonance. The results are compared in 
fig. E| with the rates given in the NACRE compilation 
|26|. They confirm the conclusion of [5| that the rate has 
to be reduced by nearly a factor of two at low tempera- 
tures, but it is in good agreement with NACRE 2(| above 
T 6 =150. 

It can be also concluded from the present analysis that 
the data above the 259 keV resonance are of crucial im- 
portance for a reliable extrapolation. This finding empha- 
sizes that one experiment alone cannot solve the problem 
at low temperatures and detailed analysis of the nuclear 
structure of 15 is required. 

A recent experimental determination of the total S- 
factor at very low energies down to E=70 keV [27j is in 
good agreement with the present R-matrix calculations. 



While extrapolation of these data to lower energies re- 
quires a detailed knowledge of the energy dependence of 
the various contributions to the total S-factor, this ex- 
periment :22] gives experimental certainty of the reaction 
rate better then 15% for Tg>90 without any extrapolation 
procedure. This clearly represents a major improvement 
in the evaluation of the reaction rate for this temperature 
regime. 

In conclusion, with the present determination of the 
reaction rates we confirm the astrophysical consequences 
in the determination of the age of the Globular Clusters 
quoted in |2J, and in the CNO solar neutrino fluxes 0G]- 
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Ex [keV] 



Branching [%] 



present work 



I2T1 



present work 



21 



5180.8 ±0.2 5183.0 ±1.0 
5240.9 ± 0.3 

6172.3 ±0.2 6173.3 ±1.7 
6791.7 ±0.2 6793.3 ±1.7 

7556.4 ± 0.6 7556.5 ± 0.4 



-yr 



17.1 ±0.2 
0.6 ±0.3 

57.8 ±0.3 

22.9 ±0.3 
1.6 ±0.1 



15.8 ±0.5 17.3 ±0.2 

57.5 ± 0.6 58.3 ±0.5 

23.2 ± 0.4 22.7 ±0.3 

3.5 ±0.6 1.70 ±0.07 



Table 1. Excitation energies of ls O states and low energy at the Fir — 259 keV resonance. 
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Fig. 14. The astrophysical Se. 17(E) factor for the capture to the 6.17 MeV state is shown from the LUNA studies (filled-in 
points), previous work [5| (open squares) and from the LENA group 17 (crosses). The insert shows the thick-target yield data 
for the 6.17 MeV secondary transition. The lines through the data points represent the results of our R- matrix fit. The dashed 
curve is scanned from |10|. 
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Table 2. S-factor for different capture transitions. S(E) in units of keV b for transitions to Ex (MeV) states. E e // denotes the 
effective center of mass energy (in units of keV) within the target £Q. 



G. Imbriani, et al: S-factor of 14 N(p,7) ls O at astrophysical energies 17 



1.1 



1.0 - 




0.5 1 1 

50 100 150 200 250 300 350 400 

T6 



Fig. 15. Reaction rate from present work is compared with that of the NACRE compilation |26|. The dashed curves represent 
the uncertainty of the present reaction rate. 



Ex [keV] 


present 
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0.25 ± 0.06 


1.55 ±0.34 


0.081^ 


0.49 ± 0.08 


0.15 ±0.07 


5183 


0.010 ±0.003 










5241 


0.070 ± 0.003 










6173 


0.08 ± 0.03 


0.14 ±0.05 


0.06l° ;S 


0.04 ±0.01 


0.13 ±0.02 


6793 


1.20 ±0.05 


1.41 ±0.02 


1.63 ±0.17 


1.15 ±0.05 


1.4 ±0.2 


total 


1.61 ±0.08 


3.20 ±0.54 


1.77 ±0.20 


1.68 ±0.09 


1.68 ±0.2 



Table 3. Values of astrophysical factor S(0) from present and previous work in units of keV b. 
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